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Fokker-Planck  Model  for  Thermal  Smoothing  in 
Inhomogeneous  Laser-Matter  Interactions 

Energy  transport  in  laser-produced  plasmas  is  an  important  issue  in  current  high 
energy  density  plasma  research,  e.g.,  inertial  confinement  fusion  (ICF)  [1,2],  fem¬ 
tosecond  laser  interactions  with  solids,  and  laboratory  simulations  of  astrophysical 
phenomena  using  lasers. 

In  ICF,  laser  energy,  absorbed  in  the  hot,  low  density  corona,  must  be  effectively 
conducted  toward  colder,  higher  density  regions  of  the  target.  A  large  inward  heat 
flux  from  the  critical  surface,  toward  higher  density  regions,  implies  that  newly 
heated  material  at  the  surface  of  the  cold  target  will  expand  outward  rapidly  lead¬ 
ing  to  large  fuel  compression.  Large  compressions  can  lead  to  higher  fusion  gains. 
A  limitation  in  ICF  is  pellet  degradation  and  disruption  due  to  the  Rayleigh- Taylor 
(RT)  instability  where  a  hot,  low-density  plasma  is  used  to  compress  and  accelerate 
a  cold,  high  density  plasma.  The  RT  instability  can  be  seeded  by  laser  nonunifor¬ 
mities  from  single  or  multiple  beams  and  target  surface  irregularities  [1,2].  Single 
beam  nonuniformities  are  the  most  detrimental  to  RT  seeding  and  contribute  to  im¬ 
print  through  perturbations  in  the  laser  driven  shock  front  and  acceleration  in  the 
post  shock  region.  Short  wavelength  laser  nonuniformities  <51/1  may  be  reduced  [1,2] 
due  to  thermal  diffusion  and  smoothing  such  that  <5P/P  <  <51/1  where  5P/P  is  the 
pressure  nonuniformity  at  the  ablation  surface  which  can  seed  the  RT  instability. 
However  longer  wavelength  laser  nonuniformities  may  not  be  reduced  by  thermal 
diffusion  alone  [3,4],  In  order  to  reduce  the  amplitudes  of  the  laser  uniformities, 
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techniques  for  laser  optical  smoothing  have  been  developed,  e.g.,  induced  spatial 
incoherence  [5],  spectral  dispersion  [6],  and  distributed  phase  plates  [7]. 

Models  for  nonuniform  laser-driven  spatial  and  temporal  heat  flow  evolution, 
over  a  large  range  of  scale  sizes,  are  important  for  ICF.  It  is  well  known  that 
that  the  measured  heat  flux  is  [8]  below  the  flux  predicted  by  classical  electron 
thermal  conduction.  Strong  inhibition  of  thermal  transport  below  classical  thermal 
conduction  seems  to  exist  in  both  long  and  short  pulse  laser  deposition  and  in 
both  high  and  low  Z  plasmas.  In  addition,  the  measured  heat  transport  seems 
to  be  reduced  both  inward  toward  higher  density  regions  and  laterally  toward  the 
exterior.  Several  models  have  been  advanced  to  explain  the  reduced  heat  flux 
in  ICF  experiments,  e.g.,  ion  acoustic  turbulence  [9],  DC  magnetic  field  effects 
[10],  and  kinetic  theory  in  steep  temperature  gradients  [11-15].  It  has  been  shown 
[16,17]  that  hydrodynamic  fluid  models  are  not  accurate  to  describe  heat  flow  when 
kAmfp  >1/80  where  27rk_1  is  the  scale  size  of  the  temperature  or  pressure  and  Amfp 
is  the  collisional  mean  free  path.  When  this  inequality  is  fulfilled,  a  kinetic  model, 
e.g.,  a  Vlasov  Fokker-Planck  (FP)  formulation,  must  be  used. 

Some  initial  work  has  been  performed  using  both  hydrodynamic  fluid  models  [3, 
19-21]  and  FP  models  [16,18],  to  study  thermal  smoothing  and  heat  flow  generated 
by  nonuniform,  inhomogeneous  laser  deposition  for  ICF  applications.  Previous  FP 
studies  [16,18]  have  been  primarily  devoted  to  the  problem  of  the  smoothing  of 


3 


electron  temperature  perturbations  with  a  single  scale  size  resulting  from  nonuni¬ 
form  laser  deposition.  The  smoothing  of  ion  density  perturbations,  resulting  from 
nonuniform  laser  depostion,  have  not  been  studied  quantitatively  using  FP  models 
over  a  range  of  scale  sizes.  These  ion  density  perturbations,  at  the  ablation  surface, 
can  seed  longer  time  scale  Rayleigh- Taylor  instabilities.  However,  laser  deposition, 
resulting  from  optically-smoothed  lasers,  can  be  structured  over  a  broad  range  of 
spatial  and  temporal  scales  [5-7].  An  outstanding  problem  is  a  kinetic  model  for 
the  generation  and  evolution  of  ion  density  perturbations  when  kAmfp  >  1/80. 

The  objective  of  this  paper  is  to  study  ion  density  perturbations,  using  a  fully 
kinetic  FP  model,  resulting  from  simulated  nonuniform  optically-smoothed  laser 
deposition  for  a  wide  range  of  scale  sizes.  It  is  found  that  significant  smoothing  of 
ion  density  perturbations  can  be  achieved  on  hydrodynamic  time  scales  for  a  range 
of  spatial  scale  sizes. 

It  is  assumed  that,  to  lowest  order,  the  electron  and  ion  distribution  functions  f 
are  weakly  anisotropic.  As  a  result,  we  make  the  standard  decomposition  [22]  with 
fa  =  foa  +  (v/v)  •  fia  with  a  =  e,  i.  For  the  electrons,  in  the  frame  of  the  ions: 


where  Yee  =  47r(e2/47rme)2lnAee,  the  laser  heating  source  term  S  can  be  written  as 
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S  =  ( Yee  / v2 )  <9/ <%  [  (nj  / 6v)  ( Z 2 In Aei  /In Aee )  ( 1  /  w  2  +  rA )  ( eEj]  /me ) 2]  <9f0e  / <9v  and  d/dt  = 
<9/dt  +  Vj  •  V.  In  addition,  zA  =  (47r0ne4Z’*/m2n3)/nA,  ^ee  =  (47rne4/m2n3)/nA, 
0=  (Z*  +  4.2/(Z*  +  0.24),  and  Z*  =  (Z2)/(Z).  Here  the  angle  brackets  denote  an 
average  over  all  ion  species.  In  Eq.  (1),  v,  is  the  average  ion  velocity  derived  from 
Eq.  (4)- (5).  The  quantities  Aie  and  A2e  have  been  previously  defined  [16,  22].  The 
ponderomotive  electric  field  Ep  =  (e/4mea;2)V|Eo|2.  In  Eq.  (l)-(2),  E  is  taken  to 
be  electrostatic  and  E0  is  the  laser  electromagnetic  electric  field.  Here  S  represents 
laser  inverse  bremstrahllung  absorption  using  the  Langdon  prescription  [23]. 

The  laser  electromagnetic  field  is  computed  using  a  wave  equation  for  the  slow 
evolution  of  the  wave  field.  The  equation  for  the  slow  Eq  can  be  written: 


d  ic2  (  d2  d2  \  1 


where  er  =  1  —  ne(z)/nc  and  e,  =  (z4;/a;o)rie(z)/nc.  In  the  transverse  x-direction 
the  wavenumbers  kx  for  the  spectral  content  of  the  laser  electromagnetic  field  are 
constructed  to  simulate  a  set  of  beamlets  generated  by  optical  smoothing.  The 
beamlets  are  then  overlapped  at  the  target  plane  using  a  lens  of  focal  length  f. 
The  diffraction  profile  at  the  target  plane  will  have  a  width  d  =  2Af/Di  where  A 
is  the  laser  wavelength  and  Di  is  the  beamlet  width  in  which  the  laser  intensity  is 
assumed  uniform.  The  transverse  wavenumbers  are  then  given  by  kxn  =  nDi/Af 
with  n=l,2, .  In  this  preliminary  study,  f=40  m,  A  =  0.25/rm,  Di  =  75  mm  so 


that  d  =  200//1H.  The  set  of  transverse  wavelengths  used  are  An  =  100/mi/n  with 


5 


n=l, . 40.  The  complex  amplitude  of  each  beamlet  is  taken  to  obey  Gaussian 

probability  statistics  with  independent  random  fluctuations  with  a  coherence  time 
rc  of  2  ps.  A  range  of  laser  rise  times  of  10  -  100  ps  were  used. 

For  the  ion  species  we  have: 


dio i  c  1  Ze  d  ,  2„  ,  n  Yu  d  (  A  „  ,  df0i\  n 

~at  +  3  V  '  f“  “  +Wm,  Tv  (”  E  '  f“)  -  ^Tv  {AuU  +  A2‘~d^)+  s‘ '  (4) 
^  +  wu  +  eJl^Im  =  tu  (5) 

at  rrij  ov 


where  Sie  =  -(me/mi)(nei/v2)9/9v((v2Te/mi)(9foj/9v  +  v3f0i) 

The  electric  field  E  in  Eq.  (1) ,  (2) ,  (4) ,  (5)  is  found  from  the  zero  current  condi¬ 
tion,  i.e.,  J  =  0.  We  have  numerically  solved  the  model  equation  set  consisting  of 
Eq.  (l)-(5).  The  2D  fully  kinetic  Fokker-Planck  code  is  solved  using  an  alternating- 
direction  implicit  method.  The  peak  vacuum  laser  intensity  is  4  x  1014W/cm2  and 
is  assumed  to  be  normally  incident  in  the  x-direction.  The  plasma  parameters  used 


are  such  that  the  ion  acoustic  transit  time  ts  ~  50  —  125ps  and  the  thermal  diffusion 


time  tth  —  20  —  lOOps. 

Fig.  1  shows  the  initial  electron  temperature  and  density  profiles  used  in  the 
simulation.  The  critical  surface  is  located  initially  at  approximately  at  z  =  100/im. 
The  ablation  surface  is  located  initially  at  approximately  at  z  =  180/jin.  The  initial 
electron  temperature  Te  is  taken  to  be  230  eV  with  the  initial  ion  temperature 
Ti  =  0.  The  initial  temperature  and  density  is  assumed  to  be  uniform  in  the 
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transverse  x-direction. 

Fig.  2  shows  the  longitudinal  and  transverse  intensity  profiles  used  to  simulate 
optically-smoothed  laser  deposition.  The  longitudinal  structure  is  consistent  with 
an  Airy  pattern  while  the  transverse  structure  is  random  both  in  space  and  time 
and  is  due  to  simulated  optical  smoothing.  These  intensity  profiles  are  used  in  the 
source  term  in  Eq.  (1).  A  single  Gaussian  beam  with  a  full  width  at  half  maximum 
of  80  /im  is  used. 

Fig.  3  shows  the  2D  variation  of  the  FP  electron  temperature  Te  at  t=420  ps 
with  a  laser  coherence  time  of  2  ps.  The  evolution  of  the  electron  temperarure  also 
shows  both  longitudinal  and  transverse  structure.  It  is  found  that  ponderomotive 
effects  are  important  only  in  the  underdense  plasma,  i.e.,  for  z  <  lOO/un.  Significant 
smoothing  of  the  electron  temperature  perturbations  is  seen  from  the  critical  to 
ablation  surface. 

Fig.  4  shows  the  2D  variation  of  the  ion  density  nt  at  t=420  ps.  It  is  found 
that  structure  at  and  near  the  critical  surface  develops  in  the  ion  density.  This 
structure  is  caused  by  a  hydrodynamic  instability  [24]  due  to  a  combination  of 
a  pressure  gradient  and  a  density  gradient  and  is  related  to  a  buoyancy-driven 
Rayleigh- Taylor  like  process.  If  the  ion  density  perturbation  is  written  as  5n,/no  oc 
exp(ik„x  +  yt)  then  a  hydrodynamic  instabilty  will  grow  with  growth  rate  7  = 
[(Cf/'cJ.)(4kjL2/l  +  4kjL2)(l/^Lj  -  1/LpL)]*/2.  Here  1/L,  =  fflnp/cte,  1/L  = 
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<91nnj/<9z,  C,s  is  the  ion  sound  speed,  and  is  the  polytropic  index.  It  is  found 
that  this  structure  is  persistent  and  is  generated  for  a  range  of  laser  rise  times 
and  coherence  times.  This  hydrodynamic  instability  is  modulated  in  the  transverse 
direction  by  the  optically  smoothed  nature  of  the  laser  deposition.  This  structure 
may  have  important  consequences  for  thermal  self- focussing  processes.  Significant 
smoothing  is  observed  also  for  the  ion  density  perturbations. 

In  addition,  the  root-mean-square  crrms  of  the  ion  density  perturbations  was 
computed  as  a  function  of  the  longitudinal  direction.  The  root-mean- square  crrms  is 
given  by: 


®rms  (t) 


(ni)  1 


\  r  r 

—  /  dx{rii  -  { Mi))' 
Lx  j 


(6) 


where  Lx  is  the  system  size  in  the  transverse  x-direction  and  the  angle  brackets 
denote  spatial  average.  The  ion  density  perturbations  maximize  near  the  critical 
surface  and  decrease  in  magnitude  into  the  overdense  region.  The  ion  density  per¬ 
turbations  at  the  ablation  surface  <rrms(z  =  zau)  <  0.01.  The  ion  density  perturba¬ 
tions  at  the  ablation  surface  arms( z  =  z m)  oc  t"  with  n  ~  0.6  —  1  for  rc  ~  2  —  8ps. 
The  ion  density  perturbations  were  also  computed  using  a  Spitzer  fluid  code  for 
two  flux  limiters,  i.e.,  f  =  0.07,  0.1.  It  was  found  that  the  FP  model  predicts  more 
smoothing  than  a  hydrodynamic  model  for  both  flux  limiters,  i.e.,  a^s  <  Lms ■ 

In  summary,  using  a  fully  kinetic  2D  Fokker-Planck  model,  the  generation  and 


evolution  of  ion  density  perturbations  from  nonuniform  laser  deposition  in  a  plasma 


slab  has  been  studied.  It  is  found  that  significant  smoothing  of  the  ion  density 
perturbations  from  nonuniform  optically-smoothed  single  beam  laser  deposition  can 
be  achieved  on  hydrodynamic  times  scales  over  a  range  of  scale  sizes.  In  addition, 
it  is  observed  that  the  Fokker-Planck  model  predicts  more  smoothing  than  the 
hydrodynamic  Spitzer  model.  It  is  found  that  ion  density  structure  can  be  caused  by 
a  hydrodynamic  instabilty  which  is  related  to  a  buoyancy-driven  Rayleigh- Taylor- 
like  process. 

We  wish  to  thank  S.  Obenschain,  D.  Colombant,  W.  Manheimer,  and  J.  Bates 
for  useful  discussions.  This  work  was  supported  by  the  U.S.  Department  of  Energy. 
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